Charge Localization in Disordered Colossal-Magnetoresistance Manganites 
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The metallic or insulating nature of the paramagnetic phase of the colossal-magnetoresistance man- 
ganites is investigated via a double exchange Hamiltonian with diagonal disorder. Mobility edge 
trajectory is determined with the transfer matrix method. Density of states calculations indicate 
that random hopping alone is not sufficient to induce Anderson localization at the Fermi level with 
20-30% doping. We argue that the metal-insulator transtion is likely due to the formation of local- 
ized polarons from nonuniform extended states as the effective band width is reduced by random 
hoppings and electron-electron interactions. 

72.15.Gd, 72.15.Rn, 71.50.+t, 71.55.+Jv, 71.38.+i 



The recent discovery of colossal magnetoresistance 
(CMR) U in Mn-oxides Lai^A^MnOa (where A can 
be Ca, Sr, or Ba) has generated extensive interests in 
these perovskites. The essential correlation between the 
magnetization and resistivity || in these compounds was 
well explained by the double exchange (DE) mechanism 
However, the precise nature of the ferromagnetic 
metal to paramagnetic insulator transition is still not 
well understood. Recently, Millis et al Q argued that 
Jahn- Teller ( JT) effects [^J have to be included to explain 
the observed resistivity behavior and the magnetic tran- 
sition temperature. Millis et al || and Bishop and co- 
workers then proposed that the metal-insulator transi- 
tion (MIT) is a consequence of the large to small polaron 
transition, induced by the reduction of effective hopping 
integrals at temperatures near T c . However, there are 
competing proposals (||| which advocate that the DE 
mechanism alone can explain the magnetic properties 
and MIT in the Mn-oxides. Notably, Varma J9| argued 
that random hopping in the paramagnetic phase due to 
the DE mechanism is sufficient to localize electrons and 
induce a MIT at T c . Although there are mounting exper- 
imental evidences of lattice effects |kJ and small polaron 
dynamics above T c |ll[] , it is still not clear at the present 
whether the small polaron formation for temperatures 
close to T c is the driving force of the MIT or it is merely 
a consequence of the MIT from the DE mechanism itself. 
To investigate the origin of the MIT, a precise calculation 
of the localization effects in the DE model is needed. In 
this Communication we study the localization properties 
in the DE model using the well-developed transfer matrix 
method @. 

The DE mechanism can be described by the following 
Hamiltonian: 
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where is the electronic spin operator and Si is the 
local spin S = 3/2 of three Mn t^g d-electrons. The 
operators c irla (c| w ) annihilate (create) a mobile e g elec- 
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tron with spin a at orbital-7 (7 £ {a, b}). — c' ila Ci ia . 
The random diagonal energy e^ 7 is included to account 
for substitutional disorders The Jij-term describes 
the Hund's rule coupling between the local spin of ti g 
electrons and the e g electrons; the V-teim describe the 
on-site inter-orbital electron-electron (e-e) interactions. 
In the Mn-oxides, Jh 3> t/S. Here we neglect the on-site 
intra-orbital e-e interactions and exchange interactions, 
which is reasonable in the large Jh limit. 

We consider only the limit SJn/t -^00. As a conse- 
quence, the electronic spin at each site is parallel to the 
local spin, for the low energy states of interest. To study 
these low energy states, we can use a projection opera- 
tor V to project out the high energy states. Then the 
Hund's rule coupling interaction can be dropped and the 
hopping term is renormalized to 

Uj = t[cos(^)cos{ 6 f) + s in(|) S m(|)e^-^], (2) 

in the classical spin limit S = 00, where {Qi,4>i) is the 
polar angles of the classical spin at site-i. Both finite 
Jh and S will decrease the localization effect, while the 
inter-orbital e-e interaction V will probably increase it. 
To simplify the problem further, we will neglect the e- 
e interaction here and assume that its only important 
contribution is to rcnormalize Hamiltonian parameters. 
The Hamiltonian ([!]) is then transformed to an effective 
one-electron Hamiltonian of the form 



H 
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The Hamiltonian (||) describes the dynamics of electrons 
in a static background of classical spins, with the hop- 
ping matrix elements dependent on the spin configura- 
tion via Eq. (||) . This approximation of dynamic disorder 
by static disorder is valid in the large S limit. Note that 
we also replace the two-orbitals by one effective orbital. 
If there is no Jahn- Teller coupling, the carrier density of 
each orbital would be approximately half of the doping 
x. However, static or dynamic JT effects will compli- 
cate this issue and the exact correspondence is not well 
understood at the metallic region at present. 

In the simplified Hamiltonian electron localiza- 
tion is due to both off-diagonal and diagonal disorder. 
The off-diagonal disorder is intrinsic to the DE model 
in the paramagnetic phase, and the resistivity due to 
this spin-disordered scattering was calculated in Born- 
approximation or memory function approximation 
j|. Localization via off-diagonal disorder has not been 
studied extensively. A few calculations |lf| have been 
done for uniform distributions of real hopping integrals t 
in one and two-dimensional systems. Also there are stud- 
ies on 2D random flux systems with uniform \t\ The 
most relevant work for the present model was a study by 
Economou and Antoniou jlTj for systems with a semi- 
circle distribution of hopping integrals. Their work was 
based on the localization criteria of Economou and Lic- 
ciardello |l8|. A recent calculation |19| based on the Zi- 
man criterion |2(| was also reported. None of these cal- 
culations produce the precise mobility edge. Moreover, 
the random Berry's phases (the phase of iy in Eq.(||)) are 
not included, and their importance is difficult to assess. 

Here we investigate the localization properties of 
Hamiltonian (|^) with the transfer matrix method |l2| . 
This technique, coupled with finite size scaling analy- 
sis, produces the most reliable information about the ex- 
tended or localized nature of the eigenstates. In this 
technique, one considers a bar of length TV and cross sec- 
tion M x M. One determines the largest localization 
length Am as iV — ► oo from the smallest Lyapunov co- 
efficient of the product of the random transfer matrix 
relevant to Eq.(||). The nature of the eigenstate can be 
determined by studying the scaling property of the local- 
ization length of finite systems. For extended (localized) 
states, Xm/M increases (decreases) with increasing M. 
At the mobility edge, which separates the extended from 
the localized states, Xm/M is independent of M and this 
behavior defines the Anderson transition. 

In the paramagnetic phase, the direction of the local 
spin Si is chosen to be uniformly distributed on a sphere: 
P((f>i) = l/(2vr); P(cos(0i)) = 1/2. Once the spin con- 
figuration {Si} is specified, the nearest neighbor hopping 
integrals can be obtained via Eq.(||). The random diago- 
nal site-energy is assumed to be distributed uniformly 
between [-W/2, W/2]. In our calculation, we have used 
systems with widths M = 4 — 14 and the length N on the 



order of N = 30,000 to minimize errors. 

We first investigate the localization effects in the ab- 
sence of any diagonal disorder, W = 0, to see whether 
random hoppings with Berry phases in the paramagnetic 
phase is alone sufficient to lead to localization of a large 
fraction of electrons, as has been argued in . Our re- 
sults are presented in Fig. 1, in which we plot the ratio of 
the calculated finite size localization length Xm with M, 
as a function of M for different energies. It is clear that 
the mobility edge E c , where Xm/M should be a constant 
independent of M, is located between 3.55 < \E c \/t < 3.6. 
To determine whether the Fermi level is below or above 
E c at certain doping levels, we need to know the den- 
sity of states (DOS) of the system, which we obtain by 
directly diagonalizing the Hamiltonian matrix for a fi- 
nite size (10x10x10) cluster. The integrated DOS and 
the DOS, averaged over many (100) configurations, are 
shown in Fig. 2. This indicates that less than 0.5% of 
states are below E c = — 3.56< (Fig. 2), far less than 
the required 20-30% for the CMR system. Therefore, 
the present calculation confirms the suggestion f| that 
purely off-diagonal disorder in the DE model is far from 
sufficient to localize electrons in the CMR materials at 
the 20-30% doping level. The inefficiency of the off- 
diagonal disorder can be understood by looking at the 
distributions of the amplitude of the hopping integral tij , 
P(\tij\) — 2\tij\/t 2 . This distribution has small weight for 
small \tij \ which are most important for localization. The 
presence of the Berry phase, on the other hand, weak- 
ens further the localization effect by breaking the time- 
reversal symmetry pl[ ] . In fact the mobility edge without 
the Berry phase is located at |i? c | w 3.3i. The inefficiency 
of the pure DE model to localize electrons clearly points 
towards the necessity to include other effects, such as 
JT electron-phonon coupling, electron-electron interac- 
tion effects, and diagonal substitutional disorders. 

To investigate the effect of diagonal disorders, we have 
calculated Am for different values of W and E, similar 
to the W=0 case, to determine the location of the mo- 
bility edge at fixed W. The mobility edge trajectories in 
the W-E plane is shown as in Fig. 3. The shape of the 
mobility edge bears remarkable resemblance to the mobil- 
ity edge trajectory in the Anderson model with diagonal 
disorder alone |2^] . The main difference is in the energy 
scale which is smaller in the present system due to the 
smaller average value of the hopping integrals. For the 
DE model, we obtain (\Uj\) — It, leading to an effective 
band edge ||| at E b = -4.0* at W=0. This reduction of 
band width also accounts for the smaller critical disorder 
W c at the band center. The outward shift of the mobility 
edge for small W is due to the increase of the effective 
band width with W. If E is normalized with the effective 
band width instead oft, the region of extended states will 
always shrink with increasing disorder W. In the same fig- 
ure, we have also plotted the equal-integrated-DOS lines. 
Figure 3 shows that to achieve localization of 20-30% of 
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the electrons, the presence of a substantial amount of di- 
agonal disorders W, in the range of 10-12 in units of t, is 
required. 

The ture mobility edge is controlled by the formation 
of localized polarons in systems with strong electron- 
phonon couplings. The presence of disorder changes the 
character of the polaron formation in three dimensions. 
In ordered systems, an abrupt change of the polaron 
state from nearly free type (large polaron) to self-trapped 
type (small polaron) occurs as the electron-phonon cou- 
pling reaches a critical value |p3 . In disordered systems, 
however, localized polaron of intermediate sizes can form 
from the nonuniform extended electronic states above the 
mobility edge even with moderate coupling |p5| . This 
transition occurs at £ « 10 A~ 2 / 3 , where A is the dimen- 
sionless electron-phonon coupling constant. The true mo- 
bility edge, obtained by using the coherent length £ from 
the finite size scaling analysis to be discussed below, is 
indicated in Fig. 3, for A = 0.03 and 1. For large val- 
ues of coupling constant, A >> 1, small polaron picture 
prevails. Figure 2 and 3 are the principal results of this 
work. 

The presence of the Berry phase breaks time-reversal 
symmetry, and hence the present model belongs to the 
unitary universality class. The critical property around 
the MIT is investigated using the one-parameter-scaling 
procedure Q. We have been successful (Fig. 4) in plac- 
ing all our data on the same universal scaling curve, 
A m /M = /(£/M), where £(W,E) is the scaling param- 
eter corresponding to the infinite size localization and 
correlation length in the localized and metallic regimes 
respectively. Moreover, our scaling data fall on the scal- 
ing function (shown as solid line in fig. 4) of the stan- 
dard Anderson model with diagonal disorder alone. The 
critical exponent v for the localization and correlation 
lengths, £~ \E-E C \- V , is found to be v = 1.0±0.2, con- 
sistent with the value for the standard Anderson model 
[ 12 1 . Our result further supports previous conclusions 
[26,^7j that the critical property of the disorder- induced 
MIT is not modified by broken time-reversal symmetry 
in three dimensional systems. 

In the Mn-oxides Lai-^A^MnOa, the potential fluctu- 
ation experienced by the e g electrons due to the La 3+ and 
A 2+ ion cores, if unscreened, would amount to U ~ 0.6 
eV. This is equivalent to W » 0.6[12x(l - x)} 1/2 eV = 
0.95 eV at x=0.3. Using Gutzwiller approximation fjj the 
band narrowing from e-e interactions at V = 20t is esti- 
mated to be around 0.5 for x — 0.2 and 0.6 for x = 0.3, if 
the two orbitals are degenerate. Hence W/t e ff can be as 
large as 13, if we assume t=0.15 eV. Screening certainly 
will reduce this value, so our estimate is only meant to 
stress the importance of disorder in MIT in Mn-oxides. 
We feel that disorder in CMR materials is not strong 
enough to induce localization by itself, but its sizable 
presence can certainly modify small polaron formation 
picture |,| for systems with strong but less than criti- 



cal value of electron-phonon coupling. This mechanism 
also has the advantage of not restricting to any particular 
form of electron-phonon coupling, and thus may apply to 
a wide range of CMR materials described by the double 
exchange mechanism. 

In summary, we have investigated the charge local- 
ization properties in the Mn-oxides in the classical spin 
limit. We find that a substantial amount of diagonal dis- 
order is required to localize 20-30% of electronic states. 
This suggests that large electron-phonon coupling and 
polaronic effect are necessary to explain the MIT close 
to T c . 

We would like to thank S.A. Trugman and H. Roder 
for helpful discussions. Ames Laboratory is operated for 
the U.S. Department of Energy by Iowa State University 
under Contract No. W-7405-Eng-82. This work was sup- 
ported by the directorate for Energy Research, Office of 
Basic Energy Sciences. Work at Los Alamos is performed 
under the auspices of the U.S. DOE. 



[1] R. M. Kusters et al, Physica (Amsterdam) 155B, 362 
(1989); K. Chahara et al, Appl. Phys. Lett. 71, 1990 
(1993); R. von Helmholt et al, Phys. Rev. Lett. 71, 2331 
(1993); S. Jin et al, Science 264, 413 (1994). 

[2] G.H. Jonker and J.H. Van Santen, Physica 16, 337 
(1950); ibid 19, 120 (1950); J. Volger, Physica 20, 49- 
66 (1954); E.D. Wollan and W.C. Koehler, Phys. Rev. 
100, 545 (1955). 

[3] C. Zener, Phys. Rev. 82, 403 (1951); P . W. Anderson 
and H. Hasegawa, Phys. Rev. 100, 675 (1955); P.G. de 
Gennes, Phys. Rev. 82, 403 (1951). 

[4] A. J. Millis, P. B. Littlewood, and B. I. Shraiman, Phys. 
Rev. Lett. 74, 5144 (1995). 

[5] J.B. Goodenough, Phys. Rev. 100, 564 (1955); G. Mat- 
sumoto, J. Phys. Soc. Japan, 29, 615, (1970). 

[6] A. J. Millis, B. I. Shraiman, and R. Mueller, Phys. Rev. 
Lett. 77, 175 (1996); and A. J. Millis et al Phys. Rev. B 
(1996) 

[7] H. Roder, Jun Zang, and A. R. Bishop, Phys. Rev. Lett. 

76, 1356 (1996); Jun Zang, A. R. Bishop, and H. Roder, 

Phys. Rev. B. 53, R8840 (1996). 
[8] N. Furukawa, J. Phys. Soc. Jpn. 63, 3214 (1994). 
[9] C. M. Varma, Phys. Rev. B 54, 7328 (1996). 
[10] TA. Tyson et al, Phys. Rev. B 53, 13985 (1996); P. Dai 

et al, Phys. Rev. B 54, R3694 (1996); S. J. L. Billinge et 

al. Phys. Rev. Lett. 77, 715 (1996); S.G. Kaplan et al, 

Phys. Rev. Lett. 77 2081 (1996). 
[11] G. M. Zhao et al, Nature 381 676 (1996); M. Jaime et 

al, Appl. Phys. Lett. March, 1996. 
[12] For a recent review, see B. Kramer and A. MacKinnon, 

Rep. Progr. Phys. 56, 1469 (1993). 
[13] J. Mazzaferro, CA. Balsero, and B. Alascio, J. Phy. 

Chem. Soli. 46, 1339 (1985); J.M.D. Coey et al, Phys. 

Rev. Lett. 75, 3910 (1995). 



3 



[14] 
[15] 

[16] 

[17] 

[18] 

[19] 

[20] 
[21] 

[22] 

[23] 
[24] 
[25] 
[26] 
[27] 



M.E. Fisher and J.S. Langer, Phys. Rev. Lett. 20, 665 
(1968). 

C. M. Soukoulis et al, Phys. Rev. B26, 1838 (1982); C. 
M. Soukoulis and E. N. Economou, Phys. Rev. B24, 5698 
(1981); Qiming Li, C. M. Soukoulis, and E. N. Economou, 
Phys. Rev. B37, 8289 (1988); M. Inui, S.A. Trugman, E. 
Abrahams, Phys. Rev. B 49, 3190 (1994). 

For a review, see B. Huckestein, Rev. Mod. Phys. 67, 357 
(1995). 

E. N. Economou and P. D. Antoniou, Sol. State Comm. 
21, 285 (1977). 

D. C. Licciardello and E. N. Economou, Phys. Rev. Bll, 
3697 (1975). 

R. Allub and B. A lacio, Sol. State Com m. 99, 613 (1996); 
and unpublished (c 



cond-mat/9608086) 



J. M. Ziman, J. Phys. C 2, 1230 (1969). 
This behavior should be distinguished from the localiza- 
tion effect due to purely random phases in otherwise pe- 
riodic systems, for example, in the random flux problem 
that has been extensively investigated recently fllq]. 
E. N. Economou. C. M. Soukoulis, and A. D. Zdetsis, 
Phys. Rev. B30, 1686 (1984); and A. D. Zdetsis, C. M. 
Soukoulis, E. N. Economou. and G. S. Crest. Phys. Rev. 
B32, 8711 (1985). 

The true band edge, according to the Lifshitz argument, 
is located at ±(6i + W/2). 

A. Sumi and Y. Toyozawa, J. Phys. Soc. Jpn. 35, 137 
(1973). 

M. H. Cohen, E. N. Economou, and C. M. Soukoulis, 
Phys. Rev. Lett. 51, 1202 (1983). 

M. Henneke, B. Kramer, and T. Phtsuki. Europhys. Lett. 
27, 389 (1994). 

E. Hofstetter and M. Schreiber, Phys. Rev. Lett. 73, 3137 
(1994). 



FIG. 4. The universal scaling function of the double ex- 
change model for different W, E, and M. At the critical point, 
Xm/M ~ 0.6. The solid line is the universal scaling curve for 
the standard Anderson model with diagonal disorder alone. 



FIG. 1. The ratio of the finite size localization length to the 
width, Xm/M, as a function of M, for the double exchange 
model in the absence of any diagonal disorders. The mobility 
edge is located around E c — — 3.56£. All states with E < E c 
are localized. 



FIG. 2. Th integrated density of states of the double ex- 
change model in the absence of any diagonal disorder, W = 
0. The density of states (in units of t -1 ) is shown in the in- 
set. The results are obtained from exact diagonalization of 
10 x 10 x 10 clusters averaged over 100 random spin configu- 
rations. 



FIG. 3. The mobility edge trajectory of the double ex- 
change model with diagonal disorder W for electron-phonon 
coupling canstant A =0, 0.03, and 1. The dashed lines show 
equal integrated density of states lines at 10%, 20%, 30%, and 
40% fractions. The density of states and the mobility edge 
trajectory are symmetric around E=0 axis. 
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